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Abstract 



The effective theory describing infinite mass particles with a given 
velocity, has a great interest in heavy flavor physics. It has the un- 
pleasant characteristic that the energy spectrum is unbounded from 
below; this fact is the source of the problems in the formulation of the 
euclidean theory. In this paper we present an analysis of the euclidean 
effective theory, that is rather complete and has positive conclusions. 
A proof of the consistency of the euclidean theory is presented and a 
technique for the evaluation of the amplitudes in perturbation theory 
is described. We compute also the one- loop renormalization constants 
of the lattice effective theory and of the heavy-heavy current that is 
needed for the determination of the Isgur-Wise function. A variety of 
effects related to the explicit breaking of the Lorentz symmetry of lat- 
tice regularization is demonstrated. The most peculiar phenomenon 
is that the heavy quark velocity receives a finite renormalization. Fi- 
nally, we compute the lattice-continuum renormalization constant of 
the Isgur-Wise current. It is needed for the conversion of the values 
of the matrix elements computed with the lattice effective theory, to 
the values in the full theory. 
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1 Introduction 



A test of the Standard Model can be realized in the quark mixing sec- 
tor by measuring the entries of the CKM matrix and checking the unitary 
relations. The extraction of the values of the matrix elements from the ex- 
perimental data is possible only if one is able to compute the effects of the 
strong interaction on the weak process. A first principle technique for the 
evaluation of the matrix elements of the weak hamiltonian between hadron 
states is lattice QCD Since present lattice cut-offs are at most 2 -=-3 GeV, 
it is not possible to simulate directly the dynamics of heavy quarks. To cir- 
cumvent this problem, effective theories have been constructed in which the 
heavy quark masses are sent to infinity p, |^ . 

A method to get a precise determination of | Vcb \ is offered by the analysis 
of the decays : 

B ^ Z^W + l + (1) 

The hadronic matrix elements for the processes (|1|) can be safely computed 
with lattice QCD by sending the beauty and the charm masses to infinity: 
Mfe, Mc ^ cxD. In this limit the matrix elements can be expressed in terms of 
a single function, the Isgur-Wise function ^ 0, 0] : 

{D,v\V,{0)\B,v') = ^MdMb {v^ + v'^) ^vv') 
{D*, V, e I V,iO) I B, v') = -z^MdMb e,,^p e^v'^^v^ ■ v') 
{D\v,e\A^{{))\B,v') = ^M^B{e^.{l + v-v')-v,v' ■e)^{vv') (2) 

where v' and v denote respectively the h and c quark 4-velocities. 
^(f ■ v') is normalized at zero recoil, ^(1) = 1. 

The computation of the Isgur-Wise function with lattice QCD requires 
the formulation of the theory of infinite mass quarks with non zero velocity 
[|] on an euchdean lattice. 

There are also other interesting applications. The effective theory allows 
the computation with lattice QCD of the production rate of heavy mesons 
in e^e~ annihilations 0: 

e+e-^DW + D^*\ + (3) 
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For center of mass energies far away from the masses of the cc or bb reso- 
nances, the heavy quarks are produced by the electromagnetic current with 
velocities that are not appreciably changed by the hadronization. The ap- 
proximation of neglecting recoil effects for the heavy quark dynamics is equiv- 
alent to the infinite mass limit. 

Studies of the lattice effective theories have already appeared in the lit- 
erature. The analytic continuation from minkowski to euclidean space has 
been treated in ref . [|1^ . Specific properties of the euclidean effective theories 



have been discussed in ref.[jTT|. These peculiarities originate from the fact 



that the energy spectrum of the effective theory is unbounded from below; 
it is the expansion of the energy-momentum relation of a heavy quark with 
momentum p = Mv + k for small k: 

P - {u- kf 



E = ^M^+p^ = Mvo + u-k + ^-^ + (4) 

where u is the kinematical velocity, u = dx/dt = v/vq. 

By removing the energy Mvq (it is a constant in a given velocity sector) 
and neglecting 1/M terms, one gets the energy- momentum relation of the 
effective theory, 

— * 

e = u ■ k, (5) 

that is unbounded from below. 

The presence of states with negative energy is then an intrinsic property of 
the effective theory and is related to the fact that one removes the energy 
Mvq associated to a non zero 3-momentum Mv. If the heavy quark picks 
up a residual momentum k with a component antiparallel to Mv, the energy 
decreases with respect to Mvq, and one is left with negative energies in eq.(||). 
Nevertheless, the theory with v is stable because it is generated with a 
Lorentz transformation of the static theory, and the latter has not negative 
energies (m = in eq.(^). The states with negative energies are simply an 
effect of the change of reference frame and do not give rise to any instability. 

The main consequences of the energy- momentum relation (|^), as shown 
in ref.[|TT|, are: 

i) The free propagator of the heavy quark in space-time can be defined only 
introducing an ultraviolet cut-off A on the residual momenta k. 
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ii) The correlators of systems composed of the effective quark and hght de- 
grees of freedom reproduce correctly the lowest order in 1/M of the original 
theory. Considering proper observables, it is also possible to take the contin- 
uum limit A ^ oo. These results are suggested by the explicit computation 
of a heavy-light correlator in the free case. 

iii) Simple Feynman rules in 4-momentum space cannot be derived. 



This paper continues the analysis started in ref.||Tl[ and is organized as 
follows. 

In sec. 2 we report a computation that includes the interaction among the 
effective and the light quarks, of the correlator considered in ref . for the 
free case. The results corroborate the conclusion of ref. |]lT|, namely that the 
effective theory has a sensible continuum limit and correctly reproduces the 
correlations of the full theory at lowest order in 1/M. 

In sec. 3 it is presented a technique for computing amplitudes of the euclidean 
effective theory in perturbation theory. Only slight modifications with re- 
spect to ordinary field theories are needed. One still has propagators and 
vertices and the only difference with respect to an ordinary theory is that it 
is not possible to integrate any more the loops over real domains in momen- 
tum space. An ie prescription is impossible, as is instead the case for the 
static theory [1^ , and one needs an additional rule for constructing contours 
for the energy integration. 

Sec. 4 deals with the lattice regularization of the effective theory. The problem 
of the doubling of the heavy quark species is discussed and the diagrams that 
are needed for the renormalization are calculated. 

In sec. 5 we discuss the renormalization of the lattice effective theory. A 
computation of the renormalization constants at full order as is presented. 
We compute also the renormalization constant of a current of the form J = 
hyVh^i, where h^i are two heavy quark fields with 4- velocity v and v' 
respectively, and F is a generic Dirac matrix. The renormalization of this 
operator is essential for converting the values of the Isgur-Wise function 
computed with lattice QCD to the values in the full theory. 

In general, the full order as renormalization on the lattice shows a vari- 
ety of effects related to the non-covariance of the regularization. Mass and 
wave function renormalization constants depend on the velocity of the heavy 
quark, contrary to the case of a covariant regularization (like dimensional 
regularization). The renormalization constant of the current J does not de- 
pend only on f ■ f ' (the only non trivial invariant given v and f'), but on the 
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components of v and v' separately. 

A very peculiar effect is also demonstrated, related to the fact that the 
effective theory contains an additional parameter with respect to the original 
theory, namely the heavy quark 4-velocity v^. There is a finite renormaliza- 
tion of the velocity (5f , that we have computed. It is absent in a covariant 
regularization. This effect does not spoil the effective theory of physical 
meaning and the original normalization of the velocity is manteined, i.e. 
'^'r — '^'b — ^1 where is the 'bare' velocity and vr is the renormalized 
one vr = vb + Sv. 

In sec. 6 we discuss the matching of the lattice effective theory with the 
original high-energy theory. 

Sec. 7 contains the conclusions of our analysis. 



2 Consistency of the theory 

In this section we present a study of the consistency of the euclidean effective 
theory. 

We consider for simplicity a simple regularization with a cut-off A on the 
spatial momenta. The free space-time propagator of a heavy quark in motion 
along the z axis is given by |1T| (see section 4 for a derivation): 



/A rj^t, _ _^ 

0(f) gA(iz-Mt) _ g-A(t2-Mt) 

= ^6ix)A5iy)A- : (6) 

Vq ztt tz — ut 

where in the last line the velocity u has been taken along the z axis. 6\ is a 
regularized delta function. 

We consider a correlator G{t, k) of a system composed of an effective and a 
light particle. Since what matters in this context is the singularity structure 
of the amplitudes, let us consider for simplicity scalar particles. In the free 
case the correlator G^^\tJ) is given by |TT| (see fig.l): 



G^^\t, I) = J d^xe-''^-'^iH{t, X I 0)AiAp{0 \ x, t) 



J {27c)^ vo2E{mv - k) 
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where A^? is the propagator of the hght particle of mass m, P is the total 
momentum of the composite system, q = P — Mqv is the momentum of the 
meson in the effective theory and I = q — mv is the residual momentum. 
At large \ k \ , the argument of the exponential becomes 

u-k + \k\ (8) 

and it is positive for m < 1. The negative energies of the effective particle 
are compensated by the positive energies of the light particle in the states 
with high virtuality. It is possible to take the continuum limit A ^ oo. The 
argument of the exponential and the prefactor in eq.(^ are the expansion 
of the full theory expressions, and the Green function G^^\t,l) correctly 
describes the internal dynamics of the composite system in lowest order in 
1/M. At large times t the correlator is dominated by the states with the 
lowest invariant mass (that eventually become the lowest bound state in the 
interacting theory), and it behaves like [|n[]: 

G(°)(t,r) ~ e-('""«+" ') * (9) 

The 2-point function G^^\t, I) describes then a system with an infinite mass, 
velocity u, and residual momentum / (the energy Mvq is removed), as the 
result of the correct coupling of a light particle with an infinite mass particle. 

We extend these results to the case of the interacting theory, considering 
the exchange of a single scalar massless particle. The correlator G(t, I) is 
given at this order by: 

GW(t,r) = Ga{t,l) + G,{t,l) + G,{t,l) (10) 

where the indices a — c refer to the fig. 2 at the end. 
By explicit computation one derives: 

Ga{t,l) = I d^kd^p{A{k,p}exp{-[u-k + E{mv + l-k)]t) + 

{k + p) + E{mv + l~k- p)]t ) + 
{k + p) + E{mv + I- k)+ \p\]t) + 
k+\p\+E{mv + t-k-p)]t) } (11) 

where p is the momentum of the scalar particle exchanged and k is the 
momentum of the heavy particle after the interaction. The functions A — D 



+ B{k,p)exp{ — [u ■ 
+ G{k,p)exp{ — [u ■ 
+ D{k,p)exp{ — [u ■ 
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are very complicated functions of the loop momenta k and p and we do not 
report their expression. They are the correct expansion of the corresponding 
functions of the full theory. 

For large loop momenta | k \ and \ p\, the arguments of the exponentials in 
the square brackets of eq.(|l5) are positive and behave like: 

u-k+\k\, u ■ {k + 'p)+ \ k + p \, 
u ■ {k + p)+ \ k \ + \ p \, u ■ k+ \ p \ + \ k — p \ (12) 

The mechanism of compensation of the energies works also in the interacting 
theory: the negative energies of the heavy 'quark' are compensated by the 
positive energies of the light 'quark' and/or the 'gluon'. It is then possible 
to take the continuum limit without encountering other divergences than 
those of a usual field theory. The functions A and B are multiplied by the 
exponentials with 2 energies, that depend respectively only on k and k + p. 
The divergences for A ^ oo originate from the integration of A over p, and 
of B over k — p. 

At large times t the correlator ([TlD behaves like the free one (P), apart 
from ultraviolet divergences in the vertex that can be factorized in the usual 
renormalization constants. 

For the amplitudes (b) and (c) the same considerations as for the ampli- 
tude (a) hold. The explicit expression of Gi, is: 

Gbit, I) = J d^kd^p{ E{k,p) exp( - [u-k+\p-k \ +E{mv + l- p)]t ) + 

+ [F{k,p)+ tG(k,p)\e^Y>{ -[u-p + E{mv + l-p)\t)} (13) 

where p is the momentum of the heavy quark before the emission of the 
scalar, and k the momentum of the heavy quark after the emission. The 
terms with F and G are associated respectively to the mass and wave function 
renormalization of the heavy quark and will be computed more easily in sec. 4. 
Gc is given by: 

G^{tJ) = J d^kd^p { M{k,p}exp{ - [u ■ k + E{mv + I- k - p}+ \ p \]t ) 

+ [N{k,p)+ tP{k,p)]exp{ -[u-k + E{mv + t-k)]t)} (14) 

where k is the momentum of the heavy quark and p is the momentum of the 
'gluon'. 
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In lattice regularization the formulas ( ]Tl| - p!^ have obvious modifications; 
the energies of the particles are replaced by the energies in lattice regular- 
ization. One can easily check that the mechanism of compensation of the 
energies is not spoiled by lattice effects. 

Even though the computation that we have presented takes into account 
the interaction only at the lowest order, we argue that the results shown have 
a general validity. Field fluctuations do not couple to the negative energies, 
that have a kinematical origin and do not point to any inconsistency. 



3 Contour representation of the amplitudes 

In this section we derive the rules for computing amplitudes of the euclidean 
effective theory in perturbation theory. A continuum regularization with a 
cut-off A on the spatial momenta is assumed for illustrative purposes; the 
variations for the lattice case are straightforward and will be discussed in 
sec. 4. 

As a first step, let us derive a contour representation of the euclidean 
effective propagator. 

The correct propagator of the heavy quark H{t,k), as a function of time t 
and spatial momentum k, has been derived in ref . ||1 If : 

iH{t,k) = ^e-"-^^* (15) 

It is forward in time, since it has to describe particle propagation only, and 
contains the correct energy- momentum relation (|^). Because of the expo- 
nential increase with time associated with negative energy states u ■ k < 0, 
the propagator (|T3p cannot be represented as the Fourier transform of a 4- 
momentum propagator. By allowing the euclidean energy of the heavy quark 
ki to be complex, one can write: 

M f dki exp(iM) f^„. 

tH{t,k) = / — — ^ (16) 

Jc Ztc iv^ki + v ■ k 

where the contour C approaches the real line for k^ ±oo, is oriented in the 
same way, and passes below the singularity of the integrand, at k^ = iu ■ k, 
for every sign of the energy. For positive energies C can be chosen as the 



8 



real axis, and in this case formula (|T6|) reduces to the Fourier transform, 
while for negative energies C has to be moved in the lower half plane. The 
representation ([TBp can also be derived by means of a Wick rotation in the 
complex plane of k^, the energy in Minkowski space. The propagator of the 
heavy quark in Minkowski space is given by: 

iH{k,,k) = — ':rT— (17) 

voko — V ■ k + le 

There is a pole in the lower half plane, at fco = u-k—ie. To preserve the causal 
structure of the theory the Wick rotation has to be done without crossing the 
pole; for positive energies the pole stays in the right quadrant, and one can 
rotate the axes as in ordinary field theories. For negative energies the pole 
stays in the left quadrant, the rotation of the real axis has to be accompanied 
by a deformation, and this produces the contour in eq. (]T6|) . 
In the static case u = 0, the pole of the integrand in eq.(|l6|) stays in the real 
axis and one can impose the correct analytic structure by an ie prescription 



[|T^. The ie of the static euclidean theory then has the meaning of a small 
mass or positive kinetic energy to ensure the decay of the correlations. 

Let us consider now an amplitude containing a heavy quark propagator, 
for example, the self-energy graph of fig. 4 In Minkowski space the amplitude 
is proportional to 

' ' (18) 



(27r)4 v{p + k) + ie k^ - X'^ + ie 
where p is the external momentum. 

In the lower half-plane of /cq there are the heavy quark and the gluon pole, at 
ko = —po + U ■ {p + k) — ie, Ex — ie, while in the upper half plane there is the 



'antigluon' pole, at ko = —Ex + ie, where E^ = y k"^ + X"^ (see fig. 3). The real 
line separates the poles of the particles from the poles of the antiparticles. 
In making the Wick rotation one has to deform the real axis in order to keep 
the same topology. The euclidean amplitude is then proportional to 



d^k 1 1 ^^g^ 



c (27r)4 + p^)+v. (p + k) k^ + A2 

where the contour C divides the fc4-plane in two regions, one containing the 
gluon pole and the heavy quark pole, the other containing the 'antigluon' 
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pole. The integration over of (|T9|) gives: 



' 1 (20) 



VoJ (27r)3 lip^ + Ex + u-{k + p)]2Ex 

where the integration extends now to the ordinary 3-momentum space. 
By performing the integration over ko in (pISD one gets: 



^ r d^k 1 ^21) 



VoJ (27r)3 l_p^ + Ex + u-{k + p)]2Ex 

The two amphtudes are the correct continuation one of the other, i.e. they 
give the same function of the external momentum if one sets p4 = ipo- 

From the above example it is easy to derive the general rule for the 
contour of integration C of the euclidean energy k^: C must divide the k^ 
plane into two connected regions, one containing only the poles of the full and 
the effective particles, the other containing the antiparticles poles. To satisfy 
this requirement the contour C has to be deformed during the integration 
over the spatial momenta. 

For positive energies of the heavy quark the contour of integration of /c4 
can be chosen as the real axis, and for negative energies the topology must 
remain the same. The rule above can also be formulated in the following 
way: one has to integrate k^ over the real axis by assuming that the poles 
of the effective particles stay always in the region they occupy for positive 
energies. 



4 Lattice regularization 

We assume the regularization of the euclidean effective theory that has been 
proposed in ref . [|l^ , that is forward in time and symmetric in space. Consid- 
ering for simplicity a motion of the heavy quark along the z axis, the action 

5 is given by: 

iS = Vo-^Kx) [ipix) - U^{x)'^{x - t)] + 

X 

- 2y [U,{x + z)^{x + z)- f/](x)^(x - z)] (22) 
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where /I is a versor in the direction /i, and U^{x) are the hnks related to the 
gauge field by: t/^(a;) = ex.^[-igA^{x - fl/2)]. 

Let us discuss the problem of the doubling of the heavy quark species. 
The energy-momentum relation of the heavy quark on the lattice is derived by 
computing the propagator as a function of time t and the residual momentum 
k: 

iHit k) = f — ^ = ^g-(t+l)ln(l+«,sinfe,) ^23) 

^' ^ Jc 27r Vo{l- e-'''^) sin k^ Vq 
One has then: 

6=1+ Uz sin kz (24) 

The energy is zero not only at /c^ = but also at A;^ = tt, implying that the 

lattice rcgularization has produced a duplication of the low energy excita- 
tions. A rcgularization that is forward in time and in space is also affected 
by the doubling problem. In this case the propagator is given by: 

^^'(^) = -n -^k,^ ^ ■ n =irT (25) 

Vo{l — e — iVz{l — e *'^^) 

and the energy-momentum relation is: 

e' = ln[l - iuzil - e-""^)] (26) 

— * 

The energy is a complex function of k and the doubling occurs when e' is 
purely imaginary, at: 

cot(A;^/2) = -Uz (27) 

We show now by a physical argument that the doubling has not any 
significant effect in the phenomenological applications of the effective theory. 
Consider a meson composed of a effective quark Q and a light antiquark q, 

— * 

with total momentum P. We assume that the doubling has been removed 
for the light quark. The effective theory deals with the residual momentum 
k of the meson: 

k = P- Mqv (28) 
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where Mq is the heavy quark mass. 
It holds: 



k = kq + ki 



(29) 



where kq and ki denote respectively the momentum of Q and of the light 
degrees of freedom. 

Since the large mass scale Mq is removed, one expects, after renormalization, 
I I to be of the order of the hadronic scale Aqcd, that is much less than 
the lattice cut-off 1/a. 

This implies that when (fcg)^ = ±7r/a and the energy of the effective quark 
is zero, the light quark momentum ki is very near to the ultraviolet cut-off 
and its kinetic energy is very large, of the order of 1/a. The configurations 
in which the heavy quark has a momentum at the edge of the Brillouin zone 
are then suppressed because of the large energy of the heavy-light system, 
as it should be. The situation is similar to that of a light meson composed 
of a Wilson fermion (r 7^ 0) and a naive fermion (r = 0). For small meson 
momenta | P |^ 1/a, the internal dynamics is described correctly, even 
though there is a duplication of the meson species. 

The doubling has a negligible effect also in the dynamics of the transition 
of a heavy meson into an heavy meson (the dynamics of the Isgur-Wise form 
factor). Due to the change of velocity of the heavy quark after W emission, 
the typical momentum transfer between the heavy quark and the light 
degrees of freedom may be greater than Aqqd- By dimensional arguments 
one expects ~ Kqcdv ■ v' . If also this scale is assumed to be much less 
than 1/a, the typical momentum exchanges lie in a region in which lattice 
effects are negligible. 

Assuming a convention for the Fourier transform according to which 
ip{x) ~ exp(iA; ■ x), one derives from the action (^) the following Feynman 
rules: 



1 



(30) 




(31) 



tad 



^ taU e-'^ 



(32) 
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2 



where k and k' denote respectively the momenta of the incoming and out- 
coming heavy quark, Vq and Vz are the interaction vertices of the heavy quark 
with a gluon with a polarization along the time or the z axis. V^"''^ are the 
vertices of emission of two gluons, for the case of the tadpole graph [k = k'). 
It is to note that the conventions for the sign of the Fourier transform and of 
the velocity are not independent, if one wants to interpret k as the residual 
momentum of the heavy quark. If one assumes a convention according to 
which il){x) ~ exp(— i/c-x), only the sign of /C4 changes in the above Feynman 



rules (i.e. the sign in front of Vz in eq.(p2D is changed). 



In usual lattice field theories every component of the loop momentum 
k^ is integrated in the interval [— vr, +7r], i.e. exp(zfc^) is integrated along a 
unitary circle. For the effective theory the integration contour of exp(i/c4) 
has to be distorted in order to keep the poles of the effective quarks always 
in the right region of the exp(z/c4)-plane, the region of the full particle poles. 
The rule for the contour in lattice regularization is analogous to the case of 
the continuum discussed in sec. 3. 

We apply now these rules to the calculation of the amphtudes that are 
needed for the renormalization of the effective theory. The infrared diver- 
gences are regulated by a fictions gluon mass A. 
The self-energy graph of fig. 4 is given by: 

2 r d'k vle-^^'P^+'^^-v!cos'ipz + kz/2) 1 
{P) 9 ^Fj ^27r)4 ^^(1 _ e-i(fc4+P4)) + sin(A;, + Pz) A{k) ^ ' 

where the integration region is the domain [— tt, +7r]'^ x C. Cp = J^tata = 
{N^-1)/2N ioian SU{N) gauge theory, and A (A;) = 2E(l-cos A;^) + (aA)2. 
Since this integral has to be computed numerically, it is convenient to reduce 
the integration region to a real domain. Making the contour integration 
analytically, one gets: 



167r2 71 J + - 1 

z{k) e-2^P4 _ vlcos'^{pz + kz/2) 
Vq(1 — z{k)e~'^'P^) + Vz sm.{kz + Pz) 



(34) 
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where A = ELi(1 - cos A;^) + AV2 and z{k) = 1 + A- ^/{lTAf^ 
The tadpole graph of fig. 5 is given by: 

m = ^ (.oe-<- - sinp.) / (35) 

In this case there is no need to integrate over k4 because the integrand does 
not contain any effective propagator and the integration region reduces to 
the ordinary one, [— tt, +7r]^. 

The vertex correction of the local heavy-heavy current J{x) = h.^,{x)rhy/{x), 
omitting the trivial spin structure (1 + v')/2 F (1 + v)/2, is given by (see 
fig.6): 

d^k 1 



5V = -g'^Cp J 



X 



(27r)4A(fc) 

VoVqC'^'^^ — Vzv'^ cos'^{kj2) 



(36) 



[fo(l — e *'^-*) + Vz sin A;2][fo(l — e ^'^■i) _|_ gj^ kz] 

where we have taken the motion of the two heavy quarks along the z axis, 
and we have set to zero the external momenta. 
Integrating over k^ one gets: 



167r2 TTVqV'o J J(i + A)2 - 1 



(1 — z{k) + Uz sin kz){'\. — z{k) + u'z sin /c^) 

5 Lattice renormalization 

In this section we describe the one-loop renormalization of the effective theory 
on the lattice. 

The self-energy Y.{k,v) of the heavy quark is given by the sum of the 
graphs considered in section 4: 

5^(fc,'y) = A{k,v) + T{k,v) (38) 

The bare propagator is given by: 

^Hik) = — — — — - (39) 

^' VQ{l-e-'^^)+VzSmkz + Mo-Y.{k,v) ^ ' 
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where we have inserted a bare mass term Mq to compute the mass renormal- 
ization condition. 

We impose on-shell renormahzation conditions. Near the mass-shell the prop- 
agator looks like: 

'^^^^ " {ivo - X)h + (t^. - Y)h + Mo - E(0) + 0(P) 



where 



Because of lattice effects (see later), the vector (X, Y) turns out to be not pro- 
portional to the euclidean velocity {wq, Vz)- This implies that mass and wave 
function renormalizations are not sufficient for a complete renormahzation 
of the effective theory. This effect can be interpreted as a renormahzation of 
the velocity. The velocity v appearing in eq.(JD) has to be identified with a 
'bare' velocity vb-, that is modified by the field fluctuations into a 'renormal- 
ized' velocity vr = vb + 5v. By comparing the bare propagator (|^) to the 
expression in terms of the renormalized parameters 

i{vr)^K + {vr),K + MR + 0{k^y ^^^^ 
and imposing the normalization of the velocity 

{vrY = {vB? = 1, (43) 

one gets, up to order first order in a^: 

6M = -^(0) (44) 
6Z = -tvoX-v,Y (45) 
6vz = —wqVzX — VqY (46) 

where 6Z = Z — 1. 

The explicit expression for the mass renormahzation 5M is: 



167r^ nvoJ J(^i ^ A)^ - 1 (l - z + UzSink, 
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The first term comes from tlie amplitude (Q) and lias a relativistic invariant 
form for small k; it is a function of the velocity because of hard gluons. The 
second term originates from the tadpole graph and is then a lattice effect. It 
is also a function of the velocity because of the explicit factor vq. 

The mass renormalization SM is a function of the velocity Uz- It is linearly 
divergent with the ultraviolet cut-off 1/a and can be written as: 

6M = '-^^ (48) 

The numerical values of x{u) are reported in table at the end. The numer- 
ical error is at most one unit in the second decimal place. For m = one 



recovers the static value already computed in ref . fl^, |14| . At /3 = 6 the mass 
renormalization is about 17% of 1/a for m = and decreases up to 9% at 
u = 0.7. 

We note that the mass renormalization 6M in the effective theory with 
w 7^ is in effect a renormalization of the residual momentum k'^ of the heavy 
quark. Indeed the heavy quark propagator can be written in the limit a ^ 
as: 

^ ^ (49) 



vk + 6M v{k- 6Mv) 



2 



where v is the euclidean 4- velocity, v = {ivo, v) and fo = + v 
One can restore the original form of the propagator, that has a pole at = 0, 
by defining a renormalized residual momentum kn by means of the relation: 

kR = k-5Mv (50) 

This effect has a very physical interpretation. The mass renormalization of 
the static quark is given by the energy of the Coulomb-like field surround- 
ing the colour charge. For v 0, the Coulomb field moves rigidly with the 
source, and carries the 3- momentum 6Mv in addition to the energy SMvq. If 
one wants to interpret k as the fraction of the heavy quark momentum that is 
changed in the collisions, and that is zero in the absence of interactions with 
other particles, it is necessary to subtract the constant contribution from 
mass renormalization. In practise, it is not necessary to make the subtrac- 
tion (pO|), because in the effective theory the energy-momentum relation is 
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linearized, and it does not matter if the expansion point is shifted by renor- 
mahzation. The only effect of 5Mv is an additional constant decay of the 
propagator with time, according to 



instead of 



iH{t,k) ~ e-"-'^* (52) 



According to eq.(|45|) the expression for the renormalization constant of 
the field 5Z is: 

_ g'^Cp 1 r I d^k 2vl z{k) + vlu-, sin k^ ^ 



^ . d^k \vl zik) - vl cos^(A;^/2)] \z{k) - cos A: j 

' ';i + yi)2_i [1- z{k) +u, sin k,]^ 



The first term in eq. ( |53| ) is infrared finite and comes from the differentiation 
of the momentum-dependent vertices. The second term is infrared divergent, 
and the singularity is isolated with the technique introduced in ref. []T3[; the 



remaining integral is evaluated numerically. Details are given in appendix A. 
One can write: 



Z{u) = 1 + ^ [ -21n(aA)2 + e{u) ] (54) 



The coefficient of the logarithmic term, i.e. the anomalous dimension of 
the heavy quark field, is independent of the velocity. It is indeed the same 
in every regularization and does not depend on the velocity in a covariant 



regularization |jT6[. The finite term e{uz) has a non trivial dependence on the 



velocity Uz, and the numerical values are reported in the table. For Uz = 



one recovers the static value already computed in ref. O, in 



The renormalization of the velocity 6vz is given, according to eq . (^q) , by: 

6vz g'^Cp 1 r f d^k 2vQZ + VzSinkz 

- ho / = : —T- + 



167r2 TT J ^/(l + Ay - 1 I - z + Uz sin k^ 



^ . d^k [v^ z{k) - vl cos^{kzl2)\ \z{k) - cos /c j 
/(l + ^)2_i \V - z^UzSvnkzf 
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The tadpole graph does not contribute to the renormahzation of the velocity, 
because the heavy quark propagator does not enter inside the loop and then 
it is not evaluated at large momenta. The first term in eq.(|5BD is a lattice 
effect while the second one has analog in the continuum and the integrand 
vanishes for small k {z{k), cos ^ 1 for A; — > 0). 

The velocity renormahzation is a finite effect, because the infrared diver- 
gences cancel between X and Y, and it can be written as: 



F 



c{u^) (56) 



The numerical values of c{uz) are reported in the table. At (3 = Q, formula 
( |56[) gives a positive renormahzation of the velocity 6vz that increases from 
10% atu = 0.1 up to 18% for u = 0.7. 

Let us discuss now a method that allows in principle a non-perturbative 
computation of the renormahzation of the velocity. We consider a specific 
example. 

The correlator Cm (^5 u) of a meson M composed of a light quark and an 
effective quark with kinematical velocity u behaves for t — > 00 like: 

GMit,u) ~ exp(-e(u)t) (57) 

where e{u) is the binding energy of a meson with velocity u in the infinite mass 
limit, e is not a physical quantity since it contains the mass renormahzation 
of the heavy quark SM{u), that is linearly divergent and has a complicate 
dependence on the velocity. 

The correlator Gnitju) of an hyperion composed of light quarks and an 
effective quark with velocity u has a time dependence analogous to that in 



eg . (|57|) , with e{u) replaced by £'(«), the hyperion binding energy. By taking 



the ratio of the 2-point functions, the mass renormahzation contribution 



6M{u) to the binding energies cancels |17] 

GHit,u) 



GM{t,u) 



exp[ - Ae{u)t ] (58) 



where Ae{u) = t'{u) — t{u) is the difference of the binding energies. Ae('u) 
is a physical quantity and, as such, satisfies the relativistic relation: 

Ae{u) = 7(u)Ae(0) (59) 
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where '~f{u) = l/y/l — v? is the Lorentz factor. 

By measuring Ae{u) and Ae(0) with numerical simulations, one can derive 
by means of eg. (pUD the renormalized velocity vr of the heavy quark. 

Let us consider now the vertex correction to the local heavy-heavy current 
J = hyVhyi. The amplitude has already been reported in sec. 4. A general 
parametrization is the following: 



2/^ 

SV = ^—^\2(vv')r(vv')\n(a\f + d(v,v')] (60) 



where r(x) = — 1 ln[x + \Jx^ — 1]. 



the logarithmic term in eq. fpOD , has already been computed in ref. it is a 
function of v ■ f the only non trivial invariant that can be constructed with 
the velocities v and v' of the heavy quarks. The finite term d is not universal 
and in lattice regularization depends separately on the components of v and 
v' . The constant d has been evaluated numerically for the case of one static 
quark, m' = 0, and one quark moving along one axis u = u^z. The numerical 
values of d{u) are reported in the table. 

The one-loop matrix element of the current J between heavy quark states 
is then given by: 

{K\J\ K,) = 1 + ^6Z{v) + ^SZ{v') + 6V{v, v') 

= l + ^[2{vv' rivv')-l)\nia\f + f{v,v')] (61) 



where we have used eqs.(|^ ^Oj) and we have defined: 

f{v,v') = ^e{v) + ^e{v') + d{v,v') (62) 

In the normalization point v ■ v' = 1 the anomalous dimension of J vanishes 
due to the conservation of the effective current related to the flavor symmetry 

i- 

For the case of an initial static quark u ' = and a final quark moving 
along the z axis u = uz, the above matrix element reduces to 

1 + ITT? [ (- 1^ - 2) In(aA)^ + f{u) ] (63) 

167r^ u 1 — u 
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The values of f{u) are reported in the table at the end. 

Let us discuss now the on-shell renormalization of the lattice effective 



theory in the real space [113 1, instead of in momentum space WM as we have 
done up to now. The two schemes differ on the lattice and the relation 
between them has been clarified in |1T7| . 

Near the mass-shell (i.e. at large times) the self-energy can be written 
as: 

J2{k) = -SM + 5Z{ vo{l - e-'''^) + V, sin k, ) + 0{k'^) (64) 

where for simplicity we have neglected the velocity renormalization 6vz that 
is not important in this context. 

The bare propagator of the heavy quark on the lattice at order as, as function 
of time t and momentum k is then given, at large times, by: 

iHit,k) = f^e'''''{tH{k^,k) + 
J Zn 

+ iH{ki,k)[ -5M + 5Z{vo{l-e-'^^)+v,sink,)]iH{k^,k)} 

= v^(^)g-(f+i)infi+«.sinAgn I -SM{t + l) 

Vo Vo (1 + Uzsmk;,) 

= Z^exp{-(t + l)ln[l + u^sinA;^ + 5M/'yo]} (65) 

where in the last line an exponentiation that is appropriate for large t has 
been done. 

In the continuum limit a — ^ 0, the propagator (|65|) reduces to: 

iH{t,k) = Z^exp[ - {t + l){6M/vo + u ■ k) ] (66) 
Vo 

The renormalization conditions in momentum space (|4^^6D imply then that 
the field renormalization constant Z is multiplied, for an evolution of time t, 
by the exponential with t + 1 instead of t. 

The evaluation of the Isgur-Wise function on the lattice requires the com- 
putation of a 3-point function G containing two heavy quark propagators. 
According to eq. (pB|) , the correlator G contains the factors exp — (t + 1) and 
exp — {f + 1) for the times t and t' of the evolution of the two heavy quarks. 
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The bare propagator of the heavy quark with renormahzation conditions 
in the real space is given, in the hmit a — > 0, by: 

Z'^ exp[-{6M'/vo + u-k)t] (67) 



Equating the expressions (^) and (^) one gets the relation between the 
renormahzation constants in the two schemes: 

Z' = , 5M' = 5M (68) 

There is a finite difference in the wave function renormahzation constants Z 
and Z' because the mass renormahzation SM is linearly divergent, i.e. SMa 
does not vanish as a — >• 0. 

The renormahzation constant of the operator J in the real space renor- 
mahzation scheme is given by: 

{K\J\K,) = 1 + hz'iv) + ^6Z'{v') + 6Viv, v') 



2/^ 

1 + T7r4[ 2(^^ ■ v' riv ■ v') - 1) HaXf + f'{v, v') ] 



= l + ^[(-lni±^-2)ln(aA)^ + /'(^)] (69) 

where in the last line the case m ' = and u = uz has been considered. The 
values of f'{u) are reported in the table. Note that in the static limit m = 0, 
/' = 0, implying that the heavy-heavy current J is exactly conserved on the 
lattice with real space renormahzation conditions. 



6 Matching 

In this section we consider the matching of the lattice effective theory with 
the full theory in the continuum [|T^]. 

For a comparison of the theoretical rate of the decays (Q) with the experi- 
mental one, it is necessary to convert the values of the form factors computed 
with the lattice effective theory, to the values in the original, 'true', theory. 
In general, it may be useful to regularize the full theory with a scheme dif- 
ferent from the lattice one (and to renormalize it). This implies that the 
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lattice low-energy effective theory and the 'target theory' differ both in the 
high energy excitations and in the regularization. The matching process is 
then conventionally divided in two steps: 

i) Matching between the bare amplitudes of the effective theory computed 
with lattice regularization, and the bare (or renormalized) amplitudes com- 
puted with a continuum regularization. 

ii) Matching between the amplitudes computed in the full and in the effective 
theory with the same regularization (or renormalization) scheme. 

It is easy to see that both steps i) and ii), i.e. the whole matching process, 
can be done in perturbation theory. 

In step i) we assume that infrared divergences are regulated in the same way. 
Then, the two regularizations differ only in the precise way in which they 
cut-off the high-momentum modes. The difference of the amplitudes in the 
two regularizations is related to hard parton effects, i.e. to partons with 
momenta of order ^qcd- Due to asymptotic freedom, this difference 

can be safely computed in perturbation theory. 

In step ii) one has an effective theory that is an expansion of the full 
theory for momenta much less than the heavy quark mass Mq. At zero 
external momenta, i.e. in the matching point, the loop amplitudes in the 
two theories differ only for virtual momenta of the order or greater than Mq. 
Since Mq ^ Aqcd also in this case the difference can be computed with 
perturbation theory. 



In this paper we consider step i); step ii) has been considered in ref.[|l9 

23- 



The matrix element in eq. (|6^) is given in the MS scheme by: 

{K I J I K') = 1 + ^(2 - - In i±^) ln(^/A)2 (70) 



The ratio of the MS matrix element divided by the lattice matrix element 
gives the factor by which one has to multiply the values of the lattice 
simulation, to get the MS values: 

Zm = l + (^[{2--ln 1±^) Hf^af - f{u)] (71) 



We have considered the real space renormalization scheme {p7\), that appears 
more natural for the lattice matrix element of the Isgur-Wise current. 
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One sees explicitly that the dependence on the gluon mass A cancels, implying 
that soft contributions cancel in the matching. 

For a numerical evaluation of one has to select a value for 0:5; at = 2 
GeV the lattice value is smaller than the continuum one by a factor 2.7. It 
is necessary to use a unique value of 0:5, otherwise the matching constant is 
no longer an infrared safe qiiantity. One has to make a guess for the higher 
orders of in Z^. By using the lattice value ioi (3 = Q and taking 11= 1/a, 
the matching constant Zm is 1 at it = and decreases with the velocity up 
to 0.95 Ski u = 0.7. With the continuum value of cts, Zm = 0.86 Ski u — 0.7. 

We note that in the general case v ^ Q and v 7^ 0, the matching constant 
Zm depends separately on v and v' , and not only on v-v' . The matrix elements 
of J computed on the lattice with different v and v' and the same v -v', must 
be multiplied with different matching constants Zm{v, v') to cancel the effects 
of the breaking of the Lorentz symmetry. 

7 Conclusions 

We have presented a fairly complete analysis of the effective theory for heavy 
quarks in euclidean space. The results are encouraging: the theory is consis- 
tent even though the energy spectrum is unbounded from below. It is stable 
and has a sensible continuum limit that reproduces the lowest order in 1/M 
of the full theory. It is possible to compute the amplitude of an arbitrary 
process of the effective theory in perturbation theory with a set of rules that 
include (in addition to the usual Feynman rules) a prescription for a contour 
integration over the energy. 

Naive lattice regularization of the continuum euclidean action leads to 
a duplication of the heavy quark species. This phenomenon has not any 
significant effect in the applications to the phenomenology, and it is then not 
necessary to add a Wilson-like term to the heavy quark action. 

The one-loop renormalization exhibits a series of phenomena that are 
consequences both of the non-covariance of the lattice regularization and of 
the fact that the effective theory selects a preferred direction in the space- 
time, namely the heavy quark 4-velocity. 

These peculiar effects of the lattice effective theory do not pose any con- 
ceptual problem for the matching with the original relativistic theory. Ef- 
fective theories describing heavy quarks with different velocities have to be 
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considered as different theories, due to the velocity super-selection rule that 
appears in the infinite mass limit @]. All these effects of the lattice regular- 
ization can be removed with velocity- dependent counter-terms. 



A Subtraction of the infrared singularities 



In this appendix we give the formulas that have been used for the subtraction 
of the infrared singularities of the loop integrals for a\ — > 0, and for the 
numerical evaluation of the integrals. 

The singularity is isolated by subtracting and adding back to the original 
integrand it's expansion for small momenta. The difference is infrared finite 
and can be computed numerically in the limit aA — 0. The term added back 
is simpler than the original integrand and can be integrated analytically. 

In the limit of small momenta | aki |<^ 1 the singular term in the wave 
function renormalization constant, eg (l53|), reduces to: 



h{k) = ^ , ^ ^ (72) 

VTWo ^k^ + {aXy [ VA;2 + (aA)2 + MA;, ]2 

where the usual phase-space factor (yf^Ci^/lGvr^ has been omitted. 
The function h{k) can be easily integrated on a 3-sphere of radius R: 

r ^ 2 1-1-?/ 

d?kh{k) = 21n(/?/aA)2 + lnl6- -In^^ (73) 
J u 1 — u 

Note that the linearized integral in eq.(^) depends on the velocity u, since 
the domain of integration is a 3-sphere, that is not 5*0(4) invariant. 
For the numerical computation, formula (^) is then replaced by: 

1 f d^k 2vl z{k) + vluz?>va.kz ^ 



TX 



+ - I S'kl 1 bo - vlco^^{Kl2)\{z{k) - n^cosfcj ^ 

/(I + ^)2 _ 1 [1 - z(/c) +M^sinfc2]2 



7r 



1 ^^CR^l^L } + 21ni?2 + lnl6--lnl±^-21n(aA)2 (74) 

^^0 \ k\{\k \ +uk,Y u l-u 
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where the phase-space factor has been removed. 

In the infrared hmit aki — ^ 0, the integral of the vertex correction (p7|) 
reduces to: 

d^A; V ■ v' 



-^Vqv'q J ^j: 2 + (^^)2 (^fc2 + (a^)2 + ^ . A^) (^^k^ + {aXy + u'-k) 

— [ d'k^ ^ (75) 

vr J k^ + {aXy Jk^ + {a\Y + u-k 



with the usual factor removed. In the last line we have taken -u ' = 0. 
The integral above is not easy to compute analytically, and one has to make 
a further subtraction. / has the same infrared singularity of the simpler 
integral 

L = — fd'k^ ^ ^ \ ^ = -ilni±^ln(i?/aA)2 (76) 



TT 



k I 2 +(aA)2 I A; I +M ■ A; u 1 - « 



The vertex correction is then computed numerically by means of the following 
formula: 

Zi [d'k{^^ '-^ + 

77 i ,^(1 + _ 1 [1 - z{k)] [1 - z{k) + u, sin k,] 

0{R- I k I) ^ 
I A? |2 (I A? I +mA;,) 

_ _1 In In + l In l±^ln(aA)2 (77) 

M 1 — M U 1 — U 

where 6{u) is a constant evaluated numerically: 

5(u) = L-I = 

dk k ^ . vT+l? + mA; , 1 + m , 

:| m = m- | (78) 



u Jo 1 + A;2 v^l + A;2 - -uA; 1 - m 
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Table 



Numerical values of the renormalization 

constants 



u 




e{u) 


e'{u) 


d{u) 


c{u) 






.0 


19.95 


24.48 


4.53 


-4.53 




19.95 


0.00 


.1 


19.89 


24.67 


4.88 


-4.58 


11.93 


20.00 


0.13 


.2 


19.71 


25.29 


5.97 


-4.74 


12.30 


20.14 


0.51 


.3 


19.37 


26.40 


7.93 


-5.05 


13.00 


20.39 


1.18 


.4 


18.78 


28.19 


10.98 


-5.56 


14.07 


20.77 


2.19 


.5 


17.75 


30.98 


15.61 


-6.46 


15.71 


21.27 


3.62 


.6 


15.81 


35.51 


22.86 


-8.22 


18.09 


21.77 


5.48 


.7 


11.17 


44.35 


36.37 


-14.4 


20.91 


19.88 


6.02 
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FIGURE CAPTIONS 



Fig.l: correlator of a system composed of an effective 'quark' Q (see text) 
and a light antiquark q in the free case. 

Fig. 2: the same correlator as in fig.l with one gluon exchange. 

Fig. 3: Wick rotation for the self-energy graph of the effective quark Q of 
fig.4. The crosses indicate the positions of the poles of the effective quark Q, 
of the gluon g and of the antigluon 'g. 

a) Case of positive energy of Q: e > 

b) Case of e < and | e |< i^A, where Ex = \jk'^-\-\^ is the gluon energy. 

c) Case of e < and | e |> E\. 

Fig.4: self-energy graph of order as of the heavy quark. 
Fig. 5: tadpole graph of order as of the heavy quark. 

Fig. 6: vertex correction of order as of the current J describing the transition 
of a heavy quark with velovity v' into a heavy quark with velocity v. 
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